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We have computed numerically the Casimir force between two identical pistons inside a very long 
f—^ ' cylinder, considering different shapes for the pistons. The pistons can be considered as quantum 

f»0 , billiards, whose spectrum determines the vacuum force. The smooth part of the spectrum fixes 

the force at short distances, and depends only on geometric quantities like the area or perimeter 
of the piston. However, correcting terms to the force, coming from the oscillating part of the 
spectrum which is related to the classical dynamics of the billiard, are qualitatively different for 
classically integrable or chaotic systems. We have performed a detailed numerical analysis of the 
corresponding Casimir force for pistons with regular and chaotic classical dynamics. For a family 
of stadium billiards, we have found that the correcting part of the Casimir force presents a sudden 
change in the transition from regular to chaotic geometries. 

I. INTRODUCTION 

a 

In the last years, there has been an increasing interest in the Casimir effect [1|. The Casimir force, which in its 

simplest version consists in the interaction between two perfectly conducting surfaces, involves the summation of the 

vacuum energy of the electromagnetic modes allowed by those surfaces, and therefore depends in a complex way on 

the spectrum associated to the considered geometry. 

^. ■ In this paper we will consider the following system: a pair of identical pistons of arbitrary shape separated by a 

distance a, contained in a very long cylinder with the same transversal shape [2]. All surfaces will be considered 

to be perfect conductors. This system has been solved, in the sense that there is an exact formula that relates the 

(renormalized) Casimir energy or force with the spectrum of the two-dimensional Laplacian operator constrained 

to the geometries that define the pistons Q. The exact formula, first derived in [3[ (see Eq. [5] below), involves a 

summation of functions of the eigenvalues of the two-dimensional Laplacian. This summation is different from the 

usual sum over modes that defines formally the vacuum energy, since the latter involves the eigenvalues of the three 

2 i dimensional problem and gives a convergent result only after a proper subtraction. The (formal) sum over modes 

is not adequate to implement direct numerical evaluations of the Casimir energy even in the simplest situations [J] , 

^ . because of its highly oscillatory behavior. However, as we will see in many specific examples, this is not the case with 

the exact formula for pistons. 

With the possibility to compute the Casimir force using the spectra of two-dimensional pistons, considered as 
quantum billiards, let us ask a natural question: is there any signature of the classical dynamics on the billiard in 
the Casimir effect? This conceptual question is based on the fact that spectra of chaotic and non-chaotic billiards are 
qualitatively different. Whereas the spectrum of a classically chaotic system shows level repulsion and its statistical 
properties are well described by Random Matrix Theory (RMT) Q, the energy levels of a regular classical system are 
Poisson distributed [6| (see [7[ for a general discussion). Moreover, the spectrum of a given geometry can be written 
as the sum of a smooth plus an oscillating part. The smooth contribution to the spectrum depends only on geometric 
quantities like the area, perimeter, and the curvature of the border. But the oscillating part depends on quantities 
related to the classical periodic orbits of the system (as action, period and stability) , which are completely different 
in the regular or chaotic case. Moreover, periodic orbits in regular systems come in continuous families, whereas in 
chaotic systems they are isolated. As already shown in Ref. [8|, Q , the smooth part of the spectrum gives the short 
distance behavior of the force between pistons, the leading term being the well known Proximity Force Approximation 
(PFA). Therefore, any signal of quantum chaos, if present, should pop out in the difference between the full Casimir 
force and its smooth part. 

The goal of this paper is to analyze the Casimir force between pistons trying to see some signature of chaos. For this 
reason we have computed the force between pistons of different shapes with regular and chaotic classical dynamics. In 
order to compute the Casimir force in our systems, it is necessary to handle summations of functions of the eigenvalues 
of the considered geometry. For very simple geometries, like a rectangle or equilateral triangle, the eigenvalues have 
explicit analytic expressions in terms of two quantum numbers, and the computation is relatively easy (although, as 
we will see, a large number of eigenvalues are necessary in order to reproduce the PFA limit). 

For circles, the eigenvalues are defined as the zeros of Bessel functions. In this case, one can follow two different 
approaches: to compute numerically the eigenvalues, and then perform the summation, or, alternatively, to perform 
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the sum over the eigenvalues by using the argument theorem. 

Besides, for rectangular and equilateral triangular billiards, the correcting part of the Casimir force, that comes 
from the oscillating part of the spectrum, can be computed using exacts results from semiclassical periodic orbit 
theory pj|. We remark that semiclassical techniques have been already used in the computation of Casimir force 
11]. For more complex geometries, semiclassical theory is only approximated and the effort to compute the periodic 
orbits is enormous. The only alternative is to find the eigenvalues associated to a given geometry using sophisticated 
methods. We will show that, with some computational effort, the exact formula is useful to evaluate the Casimir 
interaction between pistons for both integrable and non-integrable geometries. 

The paper is organized as follows. In Section II we review the main results for the force between pistons. We include 
a discussion of the exact formula, the behavior of the Casimir energy at large distance, and the relation between the 
short distance behavior and the smooth part of the spectrum. In Section III we compute the Casimir force between 
pistons with integrable classical dynamics. We will consider rectangular, triangular and circular pistons. We will use 
the large and short distance behaviors as benchmarks for the numerical calculations, and show that, up to a given 
accuracy, it is possible to find analytic expressions for the force for any distance between the pistons. We will also 
compute the correcting part of the force, and find that circular pistons have a different behavior in the short distance 
limit: the leading contribution to the correcting part of the force is finite for triangles and rectangles, and diverges 
for circles. In Section IV we present the results for the force between pistons that are billiards with chaotic classical 
dynamics. The calculation of the Casimir force in this case was done by computing the corresponding eigenvalues 
through the scaling method [121 ] . In particular, we will analyze a family of stadium billiards, and the uniformly 
hyperbolic Sinai- type billiard introduced in Ref. [131] . For each considered geometry, we will compute the full Casimir 
force, paying particular attention to its correcting part. We will show that at short distances all geometries with curved 
boundaries present a similar divergent behavior. Section V contains an analysis of the Casimir force associated to the 
oscillating part of the spectrum for both integrable and chaotic pistons. The analysis is based on the computation of 
the correcting part of the force for a family of billiards that interpolates between a regular geometry (quarter of circle) 
and chaotic billiards. We will see that, in the short distance limit, the correcting part of the force has an abrupt 
change in the transition from the quarter of circle to a stadium billiard. We present our conclusions in Section VI. 

II. CASIMIR FORCE BETWEEN PISTONS: THE EXACT FORMULA 

We consider a very long electromagnetic cylindrical cavity, with an arbitrary section (which we will assume to be 
simply connected). The cavity contains two plates (pistons) separated by a distance a. All surfaces are perfectly 
conducting. The z-direction is the axis of the cavity, and we will denote by xi the coordinates in the transverse 
sections. At classical level, the electromagnetic field admits a description in terms of independent Transverse Electric 
(TE) and Transverse Magnetic (TM) modes, which are defined with respect to z-direction. This is possible due to 
the particular geometry we are considering, that have an invariant section along the z axis. As the section is simply 
connected, there are no Transverse Electromagnetic Modes (TEM). The relevance of TEM modes in the Casimir force 
between pistons with non-simp ly connected sections has been analyzed in Ref. [141]. 

As discussed in Ref. [14J , TE and TM Casimir energies correspond to that of a set of massive scalar fields in 1 + 1 
dimensions, with masses given by the eigenvalues associated to the two-dimensional pistons, that we will denote by 
AfcN and A/cd • More explicitly, the eigenfrequencies associated to the TE and TM modes are 



^TE _ ,/^Ul? 



= V (^) 2 + A^ (n= 1,2,3,..) 
<$? = y/i'fr + ^D (n = 0,1,2, ...),, 



where the eigenvalues are defined by 



V ± i^TE,TM = — A N D </?TE,TM ■ (2) 



The eigenfunctions v?te,tm satisfy Neumann and Dirichlet boundary conditions, respectively, on the border of the 
transversal section. The pte = const eigenfunction with An = should be excluded because it does not correspond 
to a physical electromagnetic solution. 

The Casimir energy E m and force F m for a field of mass to in 1 + 1 dimensions has been computed previously by 
many authors [15| . They are given by 

£ ™( a ) = -^ L 1 ( 3 ) 



and 

F m ( fl ) = -^ = iy m ^K[(2lma) , (4) 

1=1 

where K\ is the modified Bessel function of the second kind and the prime denotes derivative with respect to the 
argument. Using these results and the analogy between the TE and TM eigenfrequencies with the eigenfrequencies 
of massive scalar fields in 1 + 1 dimensions, we can easily obtain the TE and TM contributions to the Casimir force 
between pistons in the cylindrical cavity 



F TE (a) + F TM (a) = I ]T (^ \ 2 m K[(2lX kN a) + J2 \ 2 mK[(m k ua) J . (5) 

1 = 1 \AfcN Afc D / 

This equation has been previously obtained in Ref. [3] using a different method. We stress that the formula is valid for 
a cavity of arbitrary section. In the rest of the paper we will be mainly concerned with the particular case of Dirichlet 
boundary conditions (which correspond to TM modes), because this will be enough for our purposes. Therefore we 
will omit the subindex D in the eigenvalues and the supraindex TM in the forces. 

In the next sections we will present numerical evaluations of the Casimir force between pistons using Eq. [5j In the 
cases where the eigenvalues are known analytically or numerically, the evaluation will be performed by computing 
explicitly the summations over I and Afc in Eq. [5j There are some cases where the eigenvalues associated to a given 
geometry are known implicitly through the zeros of a function /(A) = 0. In this case, the sum over eigenvalues can 
be performed using the argument theorem 

Y^X 2 k K[(2l\ k a) = 2m I dz z 2 'K[(2zla)tff , (6) 

Xk Jc f{z) 

where the integral is performed in the complex plane z along a closed contour C that encloses all zeros of f(z). 
We will now describe some general properties of the force between pistons given in Eq. [S] The function 



K[(x) = -±(K (x) + K 1 (x)) 



is a negative and increasing function for real and positive arguments, and vanishes exponentially for x — > oo. Therefore, 
the force between pistons is always attractive and decreases with distance. Moreover, at large distances, the sum is 
dominated by the term with the lowest eigenvalue Ai and 1 = 1, that is 

1 1 / A 3 \ 1/2 

F(o)«F 00 (o) = -— A?(ifo(2A 1 o) + if 1 (2A 1 o))«--(^) e - 2Aia . (7) 
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The exponential behavior is due to the finite size of the pistons, that produce a gap in the spectrum. The above 
approximation is valid as long as aAi ^> 1. 

The behavior at short distances requires a deeper analysis. As the distance between pistons decreases, more and 
more eigenvalues give a sizable contribution to the force. Taking into account the exponential behavior of the Bessel 
functions, if in the numerical computation one neglects terms that are smaller than e~ D for a given accuracy, all 
terms with 2l\^a m i n < D should be kept in the calculation, where a m i„ is the minimum distance for which the force 
is numerically computed. 

It is useful to rewrite Eq. [5] in an integral form 

F(a) = - ]T / de c K[ {2la^l)p(e) , (8) 

n i=i J 

where p(e) — ^2 k 5(e — A 2 .) = ^ is the density of states (or energy levels). The number of energy levels N below a 
given energy e can be approximated using the so called Weyl's theorem pf, which for two dimensional billiards with 
Dirichlet boundary conditions gives 

^) = (^ ~ |^ 1/2 + x) ©(e) + N(e) , (9) 



where A is the area of the piston, P its perimeter, and x is related to the shape of the boundary through 



24 ^ V «i 7T 7 12tt 

j \ / j ->7 3 

Here «i is the interior angle of each corner and k(7j) is the curvature of each smooth section jj of the border. The 
term iV(e) contains lower order contributions that vanish in the limit e — > oo. It oscillates rapidly with energy, and is 
related with the classical periodic orbits of the billiard. (For TE modes the sign of the perimeter term changes and 
the factor \ should be replaced by % — 1 in Eq. [9]). 

Introducing the first three terms of the Weyl expansion into Eq. [8l the integral and the infinite sum over I can be 
solved, giving 

F{a)& F w {a) = -— r ^A+— — -P - - — -. (11) 

87T^a 4 327ra J Ana^ 

This result has already been shown, using different methods, in Refs. j|,||, |||. The first term is the usual PFA 
proportional to the area. There are also corrections of lower order on a related to the perimeter and curvature of 
the boundary of the billiard, although the regular or chaotic classical dynamics on the surface does not enter in this 
expression. 

III. PISTONS WITH REGULAR CLASSICAL DYNAMICS 

In this Section we use the above results in order to compute numerically the Casimir force between pistons with 
rectangular, triangular and circular shapes. Thinking these shapes as billiards, they all have regular dynamics as they 
have two constants of motion (as required for two dimensional systems). We verify that the results converge to the 
expected i*V expression Eq. [IT] for short distances and to the Foo expression, Eq. [3 for large distances. 

In order to compute the Casimir force from Eq. [5J it is necessary to obtain the eigenvalues of the Laplace equation 
for each geometry. For a rectangular piston of sides L x and L yi these eigenvalues are trivially 




(12) 

with indices n,m>l for Dirichlet boundary conditions (TM modes). For Neumann boundary conditions (TE modes) 
one of the indices can be zero (but not both at the same time) . 
For an equilateral triangle of side L 
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Km = — \/n 2 + m 2 -nm (13) 

with n > 1 and m > n+ 1 for Dirichlet boundary conditions, and n > and m > n for Neumann boundary conditions, 
excluding the term (n,m) = (0,0), see [10|, LL6| f° r further details about the equilateral triangular billiard. 

On the other hand, for a circular piston of radius R, the eigenvalues of the Laplace equation are given implicitly as 
the solutions of 

Dirichlet b.c. -► J n (RKm) = 0, (14) 

Neumann b.c. — > J' n (RX nm ) — 0, (15) 

where n > is an integer which labels the order of the Bessel function, and m refers to its m-th root. For n = 0, 
the eigenvalues are non degenerated, and for n > 1 they are double degenerated. For Neumann boundary conditions 
the first trivial zero of Jq, that is Aoi = 0, is excluded. Therefore, the Casimir force can be computed using Cauchy's 
theorem as in Eq. [5] Alternatively, the eigenvalues can be computed numerically from Eq. 1141 and then use Eq. [5] to 
calculate the force. We have used both methods (see appendix [A] for details on the samples of used eigenvalues). 



A. Casimir force 

We have computed numerically the Casimir force F(a) for the three regular pistons and compared it to the predicted 
behavior for short piston separation, Fyy(a), and long separation, ^(a). The criterion to define whether a separation 
is long or short depends on the value of a compared to the characteristic lengths of each shape. Since the result for 
all shapes are qualitatively the same, we have only plotted the results of this comparison for the square with area 
A = 1 in Fig. [1] (Observe that setting the area to unity is equivalent to use units of y/area for the distance a and 
1/area for the resulting Casimir force. All figures and results in this paper are to be understood in these units.) As 
it can be seen in the short-distance region the convergence improves as we include each term of the Weyl expression 
Eq. QTJ On the other hand, in the large-distance region, we verify the convergence to the expression in Eq. [7] which 
is a function of the smallest eigenvalue. Moreover, in this region we have also plotted the next-to-leading order large 
distance asymptotic expression (including the first two eigenvalues) to show explicitly that it is possible -and simple- 
to obtain, by choosing different approximations in each region, an analytic approximation of the Casimir force which 
agrees within a given accuracy with the numerical result. Indeed, adding a few additional terms in the large distance 
expansion is possible to describe the numerical results with analytic approximations within a 1% for any value of a. 



2.0 r 




FIG. 1: (Color online) Ratio of the numeric Casimir force between squared-shape pistons to the Weyl expression (Eq. 1111) . 
Fnum/Fw, and to the large distance asymptotic expression (Eq. [7|), Fnum/Foo. For those converging to 1 at small a, from 
bottom to top corresponds to the leading (red), next-to- next-to-leading (blue) and next-to- leading order (green) in Fw (Eq.llip. 
respectively. For those converging to 1 at large a, from top to bottom corresponds to the ratio of the numeric force to the 
large-distance expression (black) Eq.[7]and to the next-to-leading order asymptotic expression (grey) -id est, including the two 
first eigenvalues. All others shapes (rectangle, triangle, circular and also the non-integrable shapes discussed in next Section) 
behave qualitatively the same. 



In order to compare the Casimir force for the integrable pistons we have plotted in Figure [5] the numeric Casimir 
force for the different geometries while keeping constant the area A = 1 for all shapes. Given the short-distance 
expression F\v(a) we find that the differences between the forces at short distances must come from the different 
perimeters of the pistons. Whereas in the large-distance region the difference must come from the first eigenvalue 
Ai (see Eq. [7|). In table of Figure [5] we give the perimeter and first eigenvalue for the geometries under study, as 
well as the geometrical \ factor. As it can be seen, there is no cross-over of the lines in Figure [2] from large to short 
distances. However this is not a general property, since it exists a spectral theory theorem that states that if region 
R contains region B! then the lowest eigenvalue of the Laplacian operator restricted to Dirichlet boundary conditions 
on R is smaller than the lowest eigenvalue of R' |17|- Therefore, it is possible to imagine pistons A and B where A is 
contained in B, but A has a larger perimeter and, therefore, the ordering in the Casimir forces for short-distances is 
the contrary to the one for large-distances. We will see an explicit example in Section IV. Another interesting prop erty 
is that, among all possible two-dimensional shapes with a given area, the circle has the lowest first eigenvalue [17j . 
Therefore, at large distances the interaction between circular pistons is stronger than for pistons of any other shape. 

We end this Section with some comments on Neumann boundary conditions (TE modes) concerning the pistons 
studied in this Section. We have computed the Casimir force in each shape for TE modes. The ordering of the results 
for small and large distances is the opposite than for Dirichlet boundary conditions (TM modes). At short distances 
this follows from the expression of Fyy : as already mentioned after equation Eq. 1111 the term proportional to the 
perimeter changes sign for Neumann boundary conditions. On the other hand, at large distances the opposite order 
comes from the relative values of the first eigenvalue of each geometry. The first eigenvalue of the billiards considered 
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FIG. 2: (Color online) Casimir force between integrable-shaped pistons with area equal to unity. The ordering for small and 
large a may be understood from the corresponding asymptotic expressions in Eq. [11] and Eq. [7J respectively, as a function 
of the ordering in the perimeter and first eigenvalue for each geometry. From top to bottom the lines correspond to a 4 : 1 
rectangle (green), an equilateral triangle (blue), a square (red) and a circle (black), respectively. The table shows the different 
parameters of each geometry that enter into Eq. [TT] and Eq. [7] 

here, for Neumann boundary conditions, is always smaller than the first one for Dirichlet boundary conditions. 
Therefore the force decreases more slowly as a — ¥ oo. 

B. Correcting part of the force 

To analyze the effects of the classical dynamics on the Casimir force we should investigate the difference 



5F(a) = F(a)-F w (a), 



(16) 



since the Weyl terms F\y only depends on geometrical quantities. The difference 8F is related to the oscillating part of 
the spectrum and will contain information of classical periodic orbits. From Eq. 111! where the last term proportional 
to x goes as a -2 , we typically expect that 



lim a 2 SF(a) 

a-i-0 



0. 



(17) 



To see if there is a qualitatively different behavior as a — > for different type of dynamics, we study numerically the 
difference of the force for regular as well as for chaotic (next Section) geometries. 

We have plotted SF for the integrable shapes studied in Figure [3] As it can be seen, the geometries consisting only 
in straight lines and angles (square, rectangle and triangle) reach a constant for short distances. For these geometries, 
semiclassical expansion can be done in an exact way, and the short distance behavior can be obtained explicitly in 
terms of periodic orbits (see Appendix [Bj . The constant can be precisely computed and as it can be seen in the 
Figure (dashed asymptotes) matches perfectly with the numerical evaluation of SF. On the other hand, the circle 
shows a divergence in SF for a — > 0. Based on the examples studied in next sections we would conjecture that this is 
a general property for billiards with curved boundaries. 

It is worth mentioning that, had we forgotten to sum a single eigenvalue in Eq. [5] of those required by the accuracy 
sought (see discussion after Eq. [7]), then the condition Eq. [T7J would be violated. Indeed, the contribution of a given 
eigenvalue Aj, to the force is 



F k (a) = -y2\ 2 k K' 1 (2l\ k a), 



(18) 



1=1 



and a simple numerical analysis shows that Fk diverges as 1/a 2 for a — > 0. Therefore the wrongly-computed SF(a) 
would diverge as 1/a 2 violating Eq. [T7J We have verified that this is not the case for all the numerical calculations. 
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FIG. 3: (Color online) Correcting part of the force (5F = F — Fw) for integrable geometries. Solid lines from top to bottom 
at a = 1 correspond to triangle (blue), square (red), circle (black) and rectangle 4:1 (green), respectively. The dashed lines are 
the asymptotes to which the polygon geometries converge according to the explicit calculations in Appendix [Bl 



IV. PISTONS WITH CHAOTIC CLASSICAL DYNAMICS 



In this Section we compute the Casimir force between pistons with non-integrable shapes. We study a family of 
stadiums billiards and a Sinai-type billiard. A stadium consist of a rectangle of sides r and £ with a quarter of circle 
of radius r, as shown in Figure |U It has been shown that the classical motion of a particle inside a stadium is fully 
chaotic [lj|. The second kind of chaotic billiard is also shown in Figure |U It has been extensively studied in Ref. pL3J j . 
and its eigenvalues have been computed [19(. 

The computation of the Casimir force requires the knowledge of the eigenvalues of the Laplace operator with a 
given boundary condition. The eigenvalues of the chaotic billiards were computed using the scaling method [12]. This 
powerful method was proposed by Saraceno and Vergini and up to date is one of the most efficients methods to solve 
the Hclmholtz eigenvalue problem in a hard- walled billiard (with Dirichlet boundary conditions). Using this method 
it is possible to find all eigenvalues and eigenfunctions in a narrow energy range by solving a generalized eigenvalue 
problem in terms of quantities over the boundary of the cavity. The gain being that in a single computational step 
a number of accurate eigenvalues are obtained in a constant proportion to the dimension of the matrices scaling as 
O(k) where k = 2tt/X is a referential wave- number. The success of the method is essentially based on the fact that 
eigenstates are quasi-orthogonal on the boundary. The scaling method was used to numerically compute extremely 
high- lying energy levels of 2D and 31? billiards [20L l2l| . 




FIG. 4: Shape of chaotic billiards pistons used in the calculations, (a) Stadiums are a family of figures formed by the union 
of a rectangle (of sides r and £) and a quarter of circle (of radius r); once the area is fixed to unity they are labeled by the 
ratio £/r. In the Figure we show a stadium with £/r = 1 (solid) and £(r — 0.2 (dashed), (b) Uniformly hyperbolic Sinai-type 
chaotic billiard, also used in this work with area fixed to unity (see [13l | for details). 



A. Casimir force 

We have computed the Casimir force between stadium pistons with Dirichlet boundary conditions for £/r = 
1, 0.705, 0.7, 0.205, 0.2, 0.005 and 0, as well as for the Sinai-type billiard, all with area A — 1 (see Appendix 
lAl for the details in the samples of eigenvalues) . 

We have verified the correct convergence of the numerical computation to the expected short and long distance 
behavior using Eq. [TTIand Eq. [71 respectively, for all billiards. 

We present in Figure [5] the numeric Casimir force for three relevant shapes. At this level there are no qualitative 
differences with the force of integrable geometries. However, in this case, we do see a cross-over of the forces between 
the Sinai-type piston and the stadiums, since the ordering in the perimeter (see table in Figure [5]) is contrary to the 
order in the first eigenvalue and, therefore, the ordering for the asymptotic small and large distance is reversed (see 
discussion in previous Section). 
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FIG. 5: (Color online) Casimir force for different shaped chaotic pistons. For small a (lower inset) from top to bottom: 
Sinai-type billiard (blue), stadium £/r = 1 (red) and stadium £/r = 0.2 (green). For large a (upper inset) from top to bottom: 
stadium £/r — 1 (red), Sinai-type billiard (blue), and stadium £/r = 0.2 (green). Due to the relationship between the perimeters 
and the first eigenvalue there is cross-over of the forces from the small to the large distance regime. At distances larger than 
those plotted in this figure there is a crossover between the green and blue curves, as expected from the fact that the lowest 
eigenvalue of the £/r — 0.2 stadium is slightly larger than the one of the Sinai-type billiard. 



B. Correcting part of the force 



In this Section we investigate the correcting part of the Casimir force that comes from the oscillating part of the 
spectrum for the non-integrable geometries under study. In all cases we found that in the limit of short distances 
lim SF(a) = ±oo. This is contrary to the the straight-line-contours integrable pistons (rectangles and equilateral 

a— V0 

triangle), but similar to the curved-line-contour (circle). Moreover, as in the circle, for these chaotic pistons with 
curved contours 5F ~ \j a as a — ¥ 0, because it is observed that lim aSF(a) = constant. These results are plotted in 

a— >Q 

Figure O In virtue of the discussion in the previous Section, we verified that a 2 SF(a) — > as a — >■ 0, and therefore 
we have not omitted any relevant eigenvalue. 



CHAOTIC TRANSITION STUDY 



In order to look for chaotic distinctive signals we have studied F(a) and 6F(a) in a transition between chaotic to 
regular geometries. These correspond to the maximal-chaotic tjr = 1 stadium to the non-chaotic £/r = stadium, 
which corresponds to the quarter of circle. In this transition, for reasons that will became clear below, we have stepped 
in the tjr = 1, 0.705, 0.7, 0.205, 0.2, 0.005 and £/r = (quarter of circle) stadiums. 

In a first stage we have computed numerically F(a) and, as expected, we have found no differences in these 
geometries that could come from a chaotic versus non-chaotic behavior. This would have been the case if the £/r = 
stadium (quarter of circle) would have had a distinctive behavior in relation to the others £/r ^ stadiums. We have 
plotted F num (a) for some relevant stadiums in Figure [7j For small a, we have that since all stadiums have the same 
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FIG. 6: (Color online) (a) Correcting part of the Casimir force for chaotic geometries, all diverge as a — > 0. (b) a8F(a) for 
these same geometries. We find that all converge to a constant. On both figures, from top to bottom the plots correspond to 
the Sinai-type billiard (blue), the stadium l/r = 0.2 (green) and the stadium l/r = 1 (red), respectively. 





area and geometrical %- factor, F num (a) is ruled by the perimeter term in Fw(a) (Eq. 1111) . It is interesting to notice 
that, for fixed area, the perimeter as a function of l/r is not a monotone function but, instead, it has a minimum at 
l/r = 1 — 7r/4 rj 0.21. This ordering in the perimeter of the stadiums is also observed in 5F(a) for large a. 

On the other hand, the study of 6F(a) has shown a qualitative separation between the chaotic stadiums and the 
non-chaotic l/r = quarter of circle. Since the expected differences should show up for small a and SF(a) diverges 
as a — > we found suitable to study this qualitative behavior in the product a SF(a) where we avoid the divergence. 
We have found that a small variation of l/r produces a small variation of aSF(a) for l/r 76 0. On the other hand, 
for l/r ~ we have found that a 5F(a) is very sensitive to l/r. To show this behavior we have plotted in Figure [8] 
the product aSF(a) for three characteristic values of l/r and a tiny variation on each of them. As it can be seen in 
the figure, the plots for l/r 76 (l/r — 0.2, 0.205, 0.7 and 0.705) show practically no difference, whereas the plots 
for l/r — and l/r = 0.005 show a large variation. In other words, the function U(£/r) = lim aSF(a) has a smooth 

a— >0 

behavior at l/r 7^ 0, and a sudden change at l/r = 0. 



r num \y>) 



perimeter 



0.08 0.10 0.1 



-600 



-800 - 




4.15 - 



4.05 



0.155 0.156 0.157 0.158 0.159 0.K 

(a) 




l/r 



FIG. 7: (a) (Color online) Numeric Casimir force for several stadiums. The lines correspond from top to bottom to stadiums 
with l/r — 1 (red), 0.7 (green), (black dashed) and 0.2 (blue), respectively, (b) Perimeter of a stadium of area equal to 1 as 
a function of l/r. Observe the correlation between both figures, since the perimeter rules the force ordering for the equal-area 
stadiums. 



VI. FINAL REMARKS 



In this work we have studied numerically the Casimir force between two identical pistons of several shapes inside a 
cylindrical cavity. We summarize our main results and conclusions. 
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FIG. 8: (Color online) aSF(a) for different values of l/r and a tiny variation on it. From top to bottom: 0.205 (blue), 0.2 (blue 
dashed), l/r — 0.705 (green), 0.7 (green dashed), 0.005 (black) and (black dashed). The plot shows the sensitivity of SF(a) 
in l/r at l/r ~ 0, where the transition from non-chaotic to chaotic occurs. 

We have presented explicit numerical calculations for the force between pistons of different shapes. The numerical 
calculations were verified using the large and small distance behaviors, which are known analytically. Moreover, we 
have shown that it is easy to obtain accurate analytic expressions for the force in the whole range of distances. 

For the simplest geometries (rectangles and triangles) , the eigenvalues are known analytically, and the evaluation of 
the force is relatively simple. For the case of the circle, the eigenvalues are known implicitly, and we used two different 
methods to compute the sum over eigenvalues: Cauchy's theorem and a numerical evaluation of the eigenvalues. 
Finally, for the more complex geometries, the only alternative is to compute the eigenvalues numerically and then 
use the exact formula for the force. The evaluation of the force at short distances involves a very large number of 
eigenvalues, and therefore a precise and efficient method to compute thousands of eigenvalues is necessary. We have 
used a particular method developed in the context of quantum billiards. At this point we would like to stress that the 
combination of this or alternative methods to compute eigenvalues |22| with Cauchy's theorem could be particularly 
useful to compute the Casimir energy for other geometries (see for instance [23j]). 

One of the original motivations of this work was to look for signals of quantum chaos in the Casimir force for the 
piston geometry. For this reason, we evaluated the part of the force originated by the oscillating part of the spectrum. 
On the one hand, we have seen a set-apart behavior depending on whether the boundary of the piston contains curved- 
lines or not. In the first case SF diverges as a — > 0, while in the second case converges. A more detailed analysis of 
the transition between chaotic and regular geometries shows a different behavior for the non-chaotic geometry in the 
part of the force coming from the oscillatory part of the spectrum, more explicitly in the finite limit U = lim a 6F(a). 

a— >Q 

We analyzed a one parameter {l/r) family of stadiums, that starts with an integrable geometry (a quarter of circle, 
l/r = 0). The finite limit is a function of this parameter, U = U{l/r) that has a smooth behavior for l/r ^ 0, and a 
jolt at l/r = 0. 

In this paper we analyzed the eventual existence of footprints of the statistical properties of the spectrum in the 
Casimir force. There could be signals of quantum chaos in the fluctuations of the Casimir force. We hope to address 
this interesting problem in a forthcoming publication. 

This work has many numerical results which have given us clues on the behavior of the Casimir Force for different 
shapes of the pistons. It would be very interesting to explore and verify analytically these results. 



Appendix A: Eigenvalues samples 



For the different shaped pistons we have used different set of eigenvalues. The following table summarizes the 
required information to reproduce our results. 
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shape 


total eigenvalues 


maximum eigenvalue 


square 


30000 


615.95 


rectangle 4:1 


30000 


616.45 


triangle 


61075 


963.35 


circle 


19984 


549.38 


1/4-circle 


16111 


451.95 


stadium tjr = 1 


17247 


467.59 


stadium tjr = 0.705 


17254 


467.7 


stadium tjr = 0.7 


21439 


521.08 


stadium tjr = 0.205 


17258 


467.67 


stadium tjr = 0.2 


18000 


477.58 


stadium £/r = 0.005 


17257 


467.63 


Sinai-type billiard 


62076 


885.47 



Appendix B: Correcting part of the force for rectangles and triangles 

As expression ([5]) for the Casimir force, the Casimir energy can be rewritten in an integral form 

E™(a) + E TE {a) = -— £ -jj dev^i (2 W3 [/to (*) + 0n(c)] , 



(Bl) 



where 



Pd,n 



( e ) = Yl S ( e ~ A M),n) = PD,N(e) + PD,N(e) 



(B2) 



is the density of states with Dirichlet (Neumann) boundary conditions, and where we have already split the density 
into a Weyl smooth part plus an oscillating part. For both rectangular and equilateral triangular billiards, there are 
explicit exact formulas for this semiclassical expansion using Poisson summation formula from explicit expressions 
(|T2"j) and (fT3")) for the eigenvalues (see [16| for details). The smooth part has the usual form for 2D billiards 



A P 

Pd,n(c) = -^ T ^—^ + Xd,n<5(e) 



(B3) 



where A = L x L y ,P = 2(L X + L y ),xr> — 1/4, Xn = —3/4 for rectangular billiard of sides L x and L y , and A — 
V3L 2 /4:, P — 3L, xd = 1/3, xn = —2/3 for equilateral triangular billiard of side L. The contribution of this smooth 
part of the Weyl series to the Casimir energy can be evaluated analytically, giving an exact result that has already 
been discussed, see Eq. [IT] 

The oscillating parts have the following exact expressions in terms of periodic orbits 



rectangle -> /od,n(e) = t~z2 Jo(L My /e) =F ^ — p Y] 
47r * — ' ZTTve ' — ' 

M v m=\ 



[L x cos(2L x my/~e) + L y cos(2L y m,^/e)] (B4) 



triangle -> pd,n(c) = 



-E 

4?r^^ 

M 



Jo(Lmy/e) T 



4vr^^ 

v m— 1 



cos((3L/2)m v / i) (B5) 



with double index M = (Mi, M2) running over integer numbers, excluding from the sum the term (Mi, M2) = (0, 0), 



which indeed had given the area term in Weyl series. The length of periodic orbits are Im = 2yJ(MiL x ) 2 + (M2L,,) 2 
for the rectangle, and Lm = \/3(M 2 + M% + M\M-i)L for the equilateral triangle. 

We see that for both shapes, and both boundary conditions, the structure of the formulae are very similar. We 
have first to evaluate the integral of terms with Jq function, that can be solved analytically 
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de^ Kl (2ia^)ML m ^) = 2a3{[2 + ( W2fl) 2 ) 2 (B6) 

Using this into (|B1[) the resulting I sum can also be performed 

El 2n 2 a 2 + 27ra sinh(7ra) cosh(7ra) — sinh (na) it 1 -2ttq\ i-d^ 

,^ (I 2 + a 2 ) 2 = 8a 4 sinh 2 (na) = 4^ " ^ + U{e > ( 7 > 

with a = i]v[/2a. We used the approximation a>l, corresponding to the limit of small plate separation a (compared 
to the typical plate's size). The final contribution of the Jo terms to the Casimir energy is 

-^uSzr-^EiUj (B8) 

\ M M M M / 

which is the same for TE and TM modes, and both shapes. As a function of plate separation a, there are a constant 
term plus a linear term on a. 

For the terms in p with cosines, we have integrals of type 

■kI 
d eK l{ 2laV-e) oob(W3 = ^ + (flm/2o)2)3/a 0*) 

with i? = 2La; or 2L y for the two terms in the rectangular billiard, and R = 3L/2 for the equilateral triangle. Now 
the sum over I can not be done analytically, but we give an approximation 

oo 

YlT, ^^ = ^?~7T^+0(e~ a ) (BIO) 

z ^ (I 2 + a 2 ) 3 / 2 a 2 2a 3 

valid for a — Rm/2a 3> 1. The sum over index m can be done in terms of Riemann Zeta function. Finally, the 
contribution of the cosine terms to the Casimir energy of TM modes (given by the Dirichlet density of states) is 

For TE modes the contribution has exactly the opposite sign. We see again a constant plus a linear term on plate 
separation a. We remark here that on semiclassical approximation (e — > oo), usually terms coming from cosine in 
Eqs. lB4l and lB5l are ignored compared to the main terms from periodic orbit, and even the Jo function is approximated 
by its asymptotic expression for big arguments. But due to the integration from low energies, we see that both Jo and 
cosine terms give contributions of the same order on a. This is the reason why a general semiclassical approximation, 
only valid for high energies, will not be enough to compute the contributions to Casimir energy. 
Collecting all terms we have the correcting part of the Casimir energy for the rectangle 

, N L X L V v^ 1 7T / 1 1 \ L X L V ^ 1 C(3) / 1 1 \ , 

and for the equilateral triangle 

^ r^ V3L 2 ^ 1 7T v^L 2 ^ 1 C(3) mi ,, 

5E DN (a) = > — =- ± ha ^- > — =- Ta ' . (B14) 

D,NV ' 32tt ^L 3 36L 8tt 2 ^ LL 9ttL 2 V ' 

M ivl M wl 
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Deriving with respect to a we finally obtain, up to exponentially small terms, the correcting part of the Casimir 
force 

r „ tmg , e -» «_, . - !g y, 4 ± i§ (z_ + zj) (B15) 



M 



J M 



which are constants. 
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